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Abstract 

We discuss the effects that the synchrotron emitting relativistic 
electrons, which are known to exist throughout many distance scales 
in active nuclei, could have on the emission line regions which 
characterize these objects. Detailed models of both the inner, 
dense, broad line region and the outer, lower density, narrow line 
region are presented, together with the first models of the optically 
emitting gas often found within extended radio lobes. We show that, 
in all cases, if the relativistic gas which produces the synchrotron 
radio emission is mixed with the emission line region gas then 
significant changes in the emission line spectrum will result. For 
small cosmic ray densities the main effect is to strengthen lines 
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which are formed in neutral regions, where photoelectric heating 

rates are lowest. As the cosmic ray; density increases, so does both 

5 

the temperature and ionization, until the cooling peak, near 10 K is 
reached and a thermal runaway occurs. 

The implications of our results for correlations between radio 
and optical properties of nuclei are discussed. We find that the 
addition of a flux of cosmic rays to a standard broad line region 
model can effectively quench Fe II emission, in agreement with Grandi 
and Osterbrock's discovery that radio loud objects tend to be weak 
Fe II emitters. The models do not reproduce their correlation 
between radio properties and Balmer decrements. Models of low 
density narrow line region gas show that relativistic particles can 
raise the temperature in the [0 III] zone by an amount sufficient to 
explain some observed 5007/4363. ratios. Although the addition of a 
flux of cosmic rays is probably not the only way to explain these 
anomalies in the emission spgctra of broad and narrow line regions, 
it does provide both a simple and natural mechanism for powering the 
optically emitting extended regions which lie in the extended radio 
lobes far from the central engine. 

The effects of the synchrotron emitting electrons on filaments in 
the Crab Nebula are discussed in an appendix, along with a comparison 
between our calculations, which employ the mean escape probability 
formalism, and recent Hubbard and Puetter models. 

Subject headings ; B1 Lacertae objects - galaxies: nuclei - nebulae: 

Crab Nebula - nebulae: general - radiation 
mechanisms - relativity 
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I . Introduction 

In general, active galactic nuclei (AGN) such as radio and Seyfert 
galaxies, as well as the more luminous quasars and qso's, are 
characterized by emission line spectra which are the composite of 
several distinct regions (see Shields and Oke 197 5, as well as the 
review by Davidson and Netzer 1979) . Photoionization models have 
been quite successful in reproducing the general features of the 
ultraviolet,' optical, and infrared spectrum (e.g., Kwan and Krolik 
1981, Pequignot 1983). However, many members of the class of AGN's 
(the radio galaxy and quasar sub-classes) are powerful synchrotron 
radio sources, and this emission is generally taken as an indication 
of the presence of particles with relativistic energies (see 
Pacholczyk 1970) . These particles could have dramatic effects on the 
cold emission line region clouds- In general, it is not known 
whether these cosmic rays actually mix with the cold gas, although it 
is known that both optical and radio emission occur on the same 
distance scales. 

It is not our intention to challenge the standard assumption that 
photoionization is an important heating and ionization process for 
emission line regions of these objects; indeed, the entire class of 
AGNs are known to be sources of relatively hard radiation (Mushotzky 
1982) and photoionization must occur. Rather, it is our goal to 
examine the effects of a modest additional heating and ionization 
agent in the form of the relativistic particles which are known to 
exist within these sources- Our motivation in this work is the 
discovery by Grandi and Osterbrock (1978) that the optical spectra of 
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radio galaxies are systematically different from those of the radio 
quiet Seyfert galaxies. Apparently, the radio properties of the 
central engine are somehow communicated to the surrounding emission 
line gas. This raises the question; what effects will the cosmic 
rays which produce the radio emission have on the emission line 
clouds? 

We know of no discussion of the effects of such a cosmic ray flux 
upon standard photoionization models of AGN emission line regions. 
Mechanical energy input in the form of shock waves (see Dopita 1977, 
Daltabuit, MacAlpine, and Cox 1978) has been considered for both 
supernova remnants and AGN's, and a model invoking particle streams 
with subsequent charge exchange neutralization (see Ptak and Stoner 
1973) has been worked out in some detail. Baldwin et al . (1977) have 
considered some properties of a model in which a flux of cosmic rays 
interact with the emission line regions, although they were primarily 
interested in setting limits on the cosmic ray flux near quasars from 
the absence of spallation products. Krolik, McKee, and Tarter (1981) 
discussed the effects of relativistic particles on a hot intercloud 
medium, and set limits to the amount of extra heating a cold cloud 
could tolerate. Balbus and McKee (1982) have discussed the stability 
of spherical clouds to disruption by energetic particles, while 
Aldrovandi and Pequignot (1973) have discussed the general case of 
cosmic ray heating. Cosmic ray heating is also thought to be 
important in heating the hot intergalactic gas found in many clusters 
(Sofia 1973, Lea and Holman 1978, Tucker and Rosner 1983). 


The next section presents some details of the calculations we have 
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performed. We use a standard photoionization code which has been 
modified to include both heating and ionization through the 
interaction of the relatively cold gas with energetic particles. In 
section III we discuss some detailed models of the various emission 
line regions of active nuclei, and show how the cosmic rays will 
affect the photoionized gas if the two mix. The last section 
presents a discussion of the implications of our results. 

A discussion of the effects of cosmic rays on a filament in the 
Crab Nebula, and a comparison between our predictions and those of 
Hubbard and Puetter (1983) for a quasar emission line region, are 
found in appendices. 


II. Model Calculations 


Our goal is to determine the effects of a flux of high energy 
particles, such as that responsible for the observed radio emission, 
on three specific models- We begin by describing the generally 
accepted model for the production of both radio continuum and optical 
emission line production in AGN, then give some details of the method 
used to model these nebulae. 


a) The Emission Regions 

The optical spectra of most AGN show evidence of at least two, and 
perhaps three, emission line regions- The evidence for this 
breakdown is summarized by Shields (1974), as well as by Osterbrock 
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(1983). Both dynamical and energy balance arguments place the dense, 

so-called 'broad line region' (BLR), at a distance of roughly a few 

light months from the central object, while a lower density 'narrow 

line region’ (NLR ) is considerably more distant (r^lO - 10 pc). 

The fact that the NLR [0 ill] lines indicate a density inconsistent 

with the presence of the familiar [0 II] and [S II] lines (see, 

Osterbrock 1981), suggests that the NLR actually has at least two 

components (cf Netzer 1982, Pequignpt 1983). In addition to these 

prominent emission regions, careful scrutiny will often reveal faint, 

very low density, emission regions within the extended radio lobes 
4 

(r ~10 pc, cf Saslaw, Tyson, and Crane, 1977, van Breugel et. al. 
1983) . 

Radio emission also occurs on all these distance scales (cf 

Osterbrock 1983, Wilson 1982, Ulvestad, Wilson, and Sramek, 1981). 

The density of the relativistic particles producing the radio 

emission can be estimated from the minimum energy argument (Burbidge, 

* _2 -3 

Jones, and O’Dell, 1974),* and these range from N 10 cm up to 
^ 10** cm ^ for the compact components near the BLR, to 'viO ^ cm ^ 
for extended radio lobes. Although it is known that the radio 
emission tends to favor an, axial symmetry (which may or may not be 
aligned with the symmetry axis of the host galaxy, Tohline and 
Osterbrock 1982) it is not. known how this axis is positioned 
relative to the emission line regions. This uncertainty is do in 
large part to a nearly complete ignorance of the geometry and 
kinematic state of the cold gas in the inner regions. 
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b) Photoionization Calculations 

Photoionization calculations are performed with a modified version 
of the program described by Ferland and Mushotzky (1982). Briefly, a 
slab of gas is divided into a large number of zones 
(typically /w 100-200), and the equations of statistical and thermal 
equilibrium are solved for each zone. The ionizing continuum is 
modified by both extinction due to the photoelectric, bremsstrahlung , 
electron scattering, and absorption line opacities, and reemission by 
the gas (treated in the outward only formalism in the first 
iteration, and by the so-called ' lambda-iteration ' method in 
subsequent iterations, see Williams 1967. Mihalas 1978). Williams 
(1967), Davidson (1972), MacAlpine (1972), and especially Davidson 
and Netzer (1979) offer further elaboration. 

The present calculations include the 'state of the art' techniques 
and processes described by Kwan and Krolik (1981) y Henry and MacAlpine 
(1982), Halpern and Grindlay (1980), and Kallman and McCray (1982). 
Line transfer is treated by the escape probability approximation as 
developed by Kwan and Krolik (1981). Hydrogen is treated with seven 
levels (2s and 2p are treated independently) and continuum, and 
escape probabilities given by; 
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These were chosen to agree with the Monte Carlo calculations of Avery 
and House (1968, see also Hummer 1962, Capriotti 1965, Hubbard and 
Puetter 1983). The Fe II ion is treated as a five level atom with 
rates and procedures described by Colin-Souf frin et al. (1980) and 
Nussbaumer, Petti, and Storey (1981). Trapping Of resonance 
transitions of Mg II and Fe II, as well as subordinate transitions of 
Fe 11^ is also treated via the escape probability formalism. The 
helium triplets are treated as a five level atom with rates and 
assumptions as in Feldman and MacAlpine (1977, with the exception 
that collisions between states which were not connected by an 
electric dipole allowed transition are given cross sections ^ 2 times 
lower, to be consistent with He spectra of planetary nebulae, see 
Kaler, Aller, Cyzack, and Epps, 1976). Rates for collisional 
excitation from the ground state to hydrogen to the first two excited 
states are taken from Aggarwal (1983). 

In addition to direct collisional excitation, we have also 
included continuum fluorescence for many lines. In the mean escape 
probability formalism, continuum pumping will proceed at a rate 

where 7^ is the continuum occupation number, A u ^ is the spontaneous 
emission rate, g 1 and g^ are the statistical weights of the lower and 
upper states, and £ is the escape probability. This approach is 
used for pumping of Fe II levels through the ultraviolet transitions 
UV1- UV2, and UV3, and the optical [Fe X] and [Fe XIV] coronal lines 

" i 

through the xuv resonance lines at 3 50A. The importance of 
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continuum pumping in the former case has been stressed by Phillips 
(1979), and for the coronal lines in Seyfert galaxies, by Osterbrock 
(1969). We find that continuum fluorescence can be a significant 
excitation process in the models we consider. 

Great care has been taken to fully incorporate the latest atomic 
cross sections and rates. The compendium by Mendoza (1983) has been 
used to update the entire collision strength and transition 
probability database. Dielectronic recombination of C, N, and 0 
through low-lying autoionizing levels has been included with the rate 
coefficients of Nussbaumer and Storey (1983). Secondary ionization 
by suprathermal electrons has been treated by combining the 
ionization efficiencies given by Shull(1979) and Bergeron and Collin- 
Souffrin (1973). . 

c) Cosmic Ray Interactions 

Synchrotron radio sources are usually modeled in terms of an 
interaction between a magnetic field and a relativistic gas with a 
typical energy per electron of a few hundred Mev (see Pacholczyk 
1970, Longair 1981). The spectral index of the radio emission for 
radio loud active galaxies is usually (V -0.7 and this suggests 
that the electrons, which make the dominant contribution to 
synchrotron emission, have a density (per unit energy interval) given 
by N^(E) /v E (Kellerman 1966). The total electron density is 
sensitive to the lower bound of the energy distribution, which is 
typically of order 10 - 100 Mev, corresponding to relativistic 
factors of ^*'1.0 - 100 (cf Lea and Holman 1978). We define a total 
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relativistic electron density as 


* 

N 



with the lower bound set to E m ^ n = 5 Mev, corresponding to a 
relativistic factor 10. This density is only weakly sensitive to 


the upper limit E max = 10 Gev because of the strong convergence of 

the electron density function We ignore any possible cosmic ray 

proton component, and present our results as a function of the total 

* 

cosmic ray electron density >N . Since the Coulomb heating rate has 


no formal dependence on the cosmic ray mass (in the relativistic 

limit), and the energy transfer rate has only a weak dependence on 

2 _ 

the relativistic factor (c*log( ), Jackson 1975), the effect of 


including a .component of relativistic protons with the same energy 
distribution as the electrons would be to double the direct heating 
rate- ( Subrelativistic protons will heat much more efficiently than 


relativistic electrons, but these particles may only affect the outer 
skin of the cloud, as discussed below* ) m 


If the cosmic rays do indeed coexist with the optical emission 
line regions, will they manage to penetrate and heat the cold gas? 
The line emitting gas in the emission line regions of an active 
galaxies ! is generally thought to be in the form of small clumps or 
filaments. These filaments are optically thick to the low energy 
ionizing continuum. As a result, the gas has a characteristic 

ionization and thermal structure; in general both the heating and 

’ - * . • ■ 

ionization rates decrease with depth into the cloud and radiative 
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transfer effects are important. 


In our discussion we assume that the cold gas is 'optically thin' 
to the energetic electrons Serious and fundamental uncertainties 
afflict detailed treatments of the penetration of energetic particles 
into cold gas, particularly if magnetic fields are present- For 
instance, if field lines are tangled near the cloud surface then 
magnetic insulation effects can impede the flow of energetic 
particles into the gas. In the simplest case penetration is impeded 
only by ionization and heating losses resulting from two body 
collisions. In this case, the ability to heat an entire cloud is 
determined by the range of a particle, or the column density of gas 
required to stop it (see Rossi 1952). Relativistic electrons have a 
range which is given within 15% by (Berger and Seltzer 1965) 


. 0.9 


R =10 
e 


25 


lOOMlJ 


particles cm 


-2 


for a gas composed of neutral hydrogen. For comparison, the column 

2 2 -2 

density of the thickest cloud we model below is ^10 cm . (The 
range of a 100 Mev electron in a fully ionized gas, in which 
bremsstrahlung and Coulomb losses are more important than ionization, 
would be some ten times smaller/ but BLR clouds typically have 
ionization fractions H°/H + /w 10.) 

The effective column density of gas seen by a cosmic ray electron 
is increased by two factors, which together can act to impede the 
flow of energetic particles into the cloud. First, the gyro radius 
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of the electron (if a magnetic field is present) is likely to be much 
smaller than the physical dimensions of the cloud. The effect is 
that the effective column density is increased by a factor of 
esc ( of ), where o< is the pitch angle of the electron. This 

effective column density is further increased by the collective 
effects discussed below. These effects increase the strength of the 
interactions between the two gasses by large and uncertain amounts, 
and in turn decreases the range of the particles by similar amounts. 
These effects depend on details of the environment near the gas 
(particularly the geometry and magnetic field) , but may decrease the 
range by factors -**10 - 100 if the relativistic and cold gasses are 
in pressure equilibrium. The clouds we discuss below remain 
optically thin to, 6 Mev eieetrons for the larger factor, although 
these details are the most uncertain part of our calculation. 
Similarly, another fundamental uncertainty is posed by magnetic 
insulation effects, which can occur if field lines are tangled near 
the surface of a cloud. 

A possible component of subrelativistic protons would probably not 

penetrate the entire cloud. For instance, a 100 Mev proton has a 

range roughly lOOx smaller than an electon of the same energy, or 
2 3 -2 

roughly 10 cm As a result, protons with considerably lower 

energies would affect only the 'skin' of the cloud. All of these 
effects introduce uncertainties which can only be resolved with a 
complete model of both the physical environment of the emission line 
region and the cosmic ray energy spectrum. For simplicity we 
consider only the relativistic electrons and assume that even the 
lowest energy ( X <v ' 10) cosmic rays are able to penetrate freely. 



-13- 


As a result the high-energy particle densities and heating rates are 
independent of depth into the cloud. Our treatment of cosmic ray 
penetration is the best that can be done with the present 
computational resources, and seems adequate for exploratory work. 


The relativistic particles both heat and ionize the gas. Our main 
concern here is for the rate with which energy is transferred to the 
cold gas (cf Lea and Holman 1978, Ginzburg and Syrovatskii 1964). In 
the H + zone of the nebula the main interaction will be with free 
electrons. Kinetic energy is passed to the cold electrons at a rate 


£ - A4/V* e.rxn 7 l<? c 


S 


(Q 


by direct Coulomb interactions (cf Jackson 1975, Cioffi 1984, Spitzer 
1962 Ginzburg and Syrovatskii 1964, Pacholczyk 1970). Here N q is 
the thermal electron density and we have integrated over the electron 
distribution given above. 

In the H° zone the main direct interaction between the two gasses 
is through ionization of the cold gas. ' For large neutral fractions 
very little of the energy of secondary electrons goes into actually 
heating the gas (see Rossi 1952, Spitzer and Tomasko 1968); 
calculations show that secondary electrons have typical energies of 
'■'■'40 ev, and that there is roughly one ionization per 15ev deposited. 
Using the rates summarized by Ginzburg and Syrovatskii (1964, the so- 
called Bethe-Bloch formula) , we find a neutral heating rate of 

^ - (\J° 3. 7 x rp ^ °q vy s'' 
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and an H° ionization rate 

p* z /.rx,o a /^A/‘c« 3 » h 

(see also Krolik, McKee, and Tarter, 1981). This ionization rate has 
been scaled through Lotz's (1967) curves to include collisional 
ionization of heavy elements in our calculations. 

Equation 6 includes only direct Coulomb interactions between the 

cosmic rays and the cold electrons, and does not include collective 

effects such as those discussed by Scott et al. (1980). As a result 

we underestimate the heating rate by factors which can be as large as 

10~*. Actually, it usually turns out that direct cosmic ray 

ionization is less important as an ionization process than the 

heating term given by equation 6. This is because the rates for 

collisional -ionization by thermal electrons are enhanced by the 

additional heat input due to the Coulomb interactions. These thermal 

ionizations occur near threshold, while cosmic ray primary 

ionizations occur far from threshold where the ionization cross 

section is small. Because the dominant effect of including 

collective effects is to increase the factor in equation 6, which in 

turn dominates the other characteristics of the model, there will be 

a homology relation between models with different collective heating 

rates and cosmic ray densities. We introduce a scale factor (3 , 

defined as the ratio of total heating to Coulomb heating, in order to 

account for this unmodelled heat input, and present our results as a 

* 

function of the product pN . (For the case studied by Scott et al., 

f =10 5 .) 
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The rates given by eqns 6-8 can be used to estimate those values 
of the electron and photon densities for which the relativistic 
particles will provide a significant source of heating and 
ionization. If every photoionization of hydrogen destroys a 1.5 ryd 
photon, the average for a typical power law continuum, then 0.5 ryd 
of heat input will precede each recombination. The ratio of particle 
to photoelectric heat input will be; 

1 = g.r»o' sx/ / pjO* C9) 
3 o,rK a % (t<l) 1\)q^ 


where 0((T e ) is the case B radiative recombination coefficient. 
Similarly, the high-energy particles provide an ionization source 
whose importance relative to photoionization is given by 


P*". /, _ <9, OS 3 p)* 

r ~ ' u 77d 



In this expression = Q( H)^ITT is the ionizing photon density 
(s cm ) , U is the ionization parameter defined as r is the 
source - cloud separation, and we assume that one neutral zone 
ionization occurs for each ionizing photon which enters the cloud. 
The last assumption is reasonable since the photons which do enter 
the neutral zone are relatively hard and capable of producing many 
secondary ionizations. 


Synchrotron emission requires the presence of a magnetic field in 
addition to the relativistic gas. If the magnetic field permeates 
the cool gas then additional cooling by optically thin cyclotron 
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emission can be significant under some circumstances. This process 
is related to the total cooling rate by (see Fabian, Pringle, and 
Rees 1976) 


f\ C y] c./o __ 

^ *h i + 


Ve <■' 7b l * c 

O.S- hv> 0 (T e ) /v e A/* 


0 , 2 * BV e 

K)q_ 


0» 


This cooling mechanism is negligible for parameters of interest here 
(B/*10 ^ G, T~10^K, 10 10 cm ^ for the BLR). We also assume that 
the prevalent magnetic field is not significantly enhanced within the 
gas cloud. Were such the case, the gyro radius of relativistic 
particles could be small enough to inhibit their flow into the cloud 
(i.e., the magnetic insulation effects alluded to above). 


Finally, we note that, for comparison with radio observations 

which estimate the pressure due to relativistic electrons, the 

* * 5 

relationship between N and this pressure is N = 1.0x10 p re ’» where 

we use the parameters given above, and assume that the pressure is 

entirely due to fully relativistic electrons in equipartition with 

the magnetic field. We also note a constraint, studied by Balbus and 

McKee (1982), of the stability of these clouds when p re ^ approaches 

the cold gas pressure. Such clouds tend to be unstable to disruption 

by pressure gradients when they are supported entirely by non-thermal 

energetic particles. 
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III. Specific Models 
a) Broad Line Region Clouds 

Baldwin's (1977) discovery of the so-called L<x /HjB problem in the 
spectra of quasars has rejuvenated interest in photoionization models 
of these objects (see Canfield and Puetter 1980, Kwan and Krolik 
1981, Weisheit, Shields, and Tarter 1981). It has long been known 
that the energy density of the relativistic particles within compact 
radio components is comparable to that of the BLR gas (Mathews 
1974). Miley and Miller (1979) have established correlations between 
compact radio structure and Fe II emission, while Grandi and 
Osterbrock (1978) show that radio galaxies have steeper Balmer 
decrements and stronger Fe II emission than Seyfert galaxies. 
Considering a mixture of radio loud and quiet objects, Clavel and 
Joly (1983) found that Fe II emitters tended to have larger 
C III] 1909/C IV 1549 intensity ratios. Apparently, some aspects of 
the emission line formation mechanism are somehow coupled to the 
radio properties of the nucleus. 

Typically the compact radio lobes of radio loud objects have 

0 — 3 ★ 8*5 — 3 

relativistic particle densities of 10 cm < N < 10 cm (see 
Burbidge, Jones, and O'Dell 1974). Although these compact radio 
components occur at roughly the BLR radius, it is an open question 
whether the two gasses actually mix. The only previous constraint on 
the cosmic ray flux deduced from the emission line spectrum was that 


of Baldwin 

et al. (1977). 

They 

set limits 

of 

* 6-3 

N <10 cm by 

considering 

the radiation 

damage 

suffered 

by 

each cloud and 
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estimating the total exposure time. 

We reconsider the 'standard' model of Kwan and Krolik (1981) in 

detail here. As is usually the case for clouds ionized by a 

continuum with a significant x-ray component, this model is 

effectively matter bounded. The original Kwan and Krolik model had a 

23 — 3 

column density of 10 cm , which was set by integrating until 

C II] 2326 reached its observed strength. This corresponded to a 

12 

total linear thickness of 6.3x10 cm. The majority of this gas is 

fairly neutral; the hydrogen half-neutral point occurs at of depth of 
12 

only 1.0x10 cm. It is the line transfer effects which occur within 
the extended neutral zone which produces the anomalous m /H|3 ratio. 

We make the same assumptions as Kwan and Krolik (1981) in setting 

parameters of the model. The ionizing continuum is their broken 

power law, and the cloud has constant gas pressure and an inner 

hydrogen density of N = 10 cm The chemical composition is 

ti 

nearly solar, with the exception of carbon, which is depleted by a 

factor of two. We use both their ionization parameter and stopping 

criteria (e.g., radial integrations cease when the 

C II] 2326/C IV 1549 intensity ratio reaches 0.1). Our predictions 

are in good general agreement with those of Kwan and Krolik (for 
* 

small N ) despite the many changes in the atomic database which have 
occurred since their calculations were performed. 

Figure 1 shows some results of these calculations. These can be 
understood as the twofold consequences of increasing the heating and 
ionization rates through the addition of a component of high energy 
particles. In any steady state model the total intensity of all 
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emission lines produced at a particular point is set by the heating 

rate at that point. In the Kwan and Krolik model (see Fig 2), a 

large portion of the total power radiated originates in the small H + 

'skin' on the inner boundary of the filament. This follows because 

the heating rate is largest where the continuum is relatively 

★ 

unextinguished. Under normal circumstances (e.g. f negligible N ), a 
very large extended neutral zone is required to produce appreciable 
emission in lines from Mg + , Fe + , C + , and H° because of the lower 
heating rate where the radiation field has been extinguished. 

When the cosmic rays are first added they affect this outer 
neutral zone more than inner regions since this is where the original 
photoelectric heating and ionization rates were smallest. A 
consequence is that the heating rate and line emissivities in the 
'back' of the cloud increases, and the column density needed to 
produce the desired C Il]/C IV ratio decreases. (Had we maintained 
constant column density rather than constant C Il“|/C IV line ratio, 
many lines originating in the neutral zone, particularly C II] 2326, 
would achieve prodigious strengt-h. ) The Lt( /H(i ratio grows larger 
because of additional LQ( production, largely because the gas is 
growing hotter. The optical depths in nearly all lines grows 
smaller, and the absolute strength of h( 3 changes in a complicated 
manner because of both the decreasing optical depth and increasing 
temperature (see also Fig. 3 of Netzer 1975). 

★ 

As jiN increases still further the ionization and temperature 
continues to increase, and even some lines which are formed mainly in 
the H + skin (e.g., C IV 1549) strengthen relative to L<X • Some 
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lines weaken (e.g., the Fe II lines) because they are no longer 

formed by dominant stages of ionization. Eventually the additional 

heating and ionization makes the dominant contribution to the thermal 

and ionization equilibrium of the cloud and a 'nebular' solution is 

no longer found. As is typical of ionized gas, the cooling curve for 

these clouds reaches a peak in the neighborhood of 10^ K (see 

Dalgarno and McCray 1972). When the heating rate exceeds this 

maximum cooling, typically at a ratio of relativistic to cold 

electrons of /N = 10 , the gas undergoes a thermal runaway and 

reaches relativistic temperatures. (This result is in good agreement 

with Krolik et al.'s (1981) estimate of the largest additional heat 

input consistent with a nebular solution. ) This runaway actually 

* 5,5 _3 

sets a more stringent limit to N ( < p 10 cm ) than does the 
radiation damage argument of Baldwin et al . (1977). 

To consider the implications of these results for correlations 

between emission line spectra and radio properties, we give results 

of two models in detail. In addition to the model with the original 

Kwan and Krolik (1981) parameters, Figure 2 shows a model with 
*4*5-3 

pN =10 cm , and Table 1 lists detailed predictions. In 
general, many of the optically thin intercombination lines strengthen 
relative to more nearly thermalized lines (such as L°() in response 

to the additional heat input. As was mentioned above, the Fe II 

* 

emission has weakened considerably because very little Fe is present 
in the model- This is mainly because the model with the cosmic rays 
was able to reproduce the desired C II] 2326 intensity with a 
relatively small column density, so that most iron is still doubly 
ionized at the back face of the cloud (the ionization potential of C + 
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is A/ 50% larger than that for Fe + ) . The calculations do not 

reproduce the steeper Balraer decrements observed in radio loud 
objects, which Osterbrock (1983) notes could be evidence for embedded 
dust, nor do they predict the larger C IIl]/C IV ratio in Fe II 
emitters . 

An interesting prediction is that the L°< line should become more 

★ 

nearly symmetric as N increases. One of the most promising models 
for the origin of BLR clouds is the radiatively accelerated 
outflowing wind as developed by Blumenthal and Mathews (1975, 1979). 
The greatest shortcoming of this model is its prediction of an 
extremely asymmetric Lb( profile. This follows because LC< is almost 
entirely formed in the H + skin on the inner side of the BLR cloud, 

i 

and hence is preferentially radiated in the direction of the central 

object (Ferland, Netzer, and Shields 1979). This is in sharp 

contradiction with the symmetric L0( usually observed (cf, Baldwin 

and Netzer 1978). As Figure 1 shows, significant emission in the 

outward direction can be produced, largely because of enhanced LP( 

emission from the back of the cloud. This happy circumstance only 

* 5-3 

occurs when {JN is larger ( > 10 cm ) than would actually be 

allowed, considering other line ratios such as L# /Hp , however. 
(One other result of the fact that all hydrogen lines are emitted 
from both sides of the cloud is that convergence of the entire model 
was slow. Typically 4 to 6 iterations were needed.) 

Interactions between the radio components and BLR clouds could 
also explain the absence of emission lines in BL Lac objects. Models 
proposed so far to account for this lack include the absence of a 
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sufficient amount of gas to produce the lines (Shields 1978), an 

alignment with the relativistic jets pointed towards the observer 

(Blandford and Rees 1978), and an intrinsic inability for emission 

line clouds to form (Guilbert, Fabian, and McCray 1983). We point 

out another possibility; namely that relativistic particles are mixed 

with the BLR gas and cause it to emit a negligible amount of optical 

* -4*4 

radiation. Figure 1 shows that if /N > 10 the gas will 

undergo a thermal runaway. Such gas would produce little optical 
emission . 


b) Narrow Line Region 

Osterbrock (1982) speculates that the relativistic plasma within 
radio jets may explain the high temperatures often deduced from the 
[O III] 5007/4363 ratio. Heckman and Balick (1979) find temperatures 

4 

in excess of 2x10 K for NLR's m their sample of broad lxne objects, 

while Cohen and Osterbrock (1981) find even higher temperatures for 

some narrow line radio galaxies. These temperatures are surprising 

4 

since photoxonxzatxon generally produces temperatures close to 10 K 
(see Osterbrock 1974). 

We model a typical NLR as a set of clouds of constant density 
3 —3 

(N = 10 cm ) and an ionization parameter (defined as the ratio of 

ionizing photons to free electrons at the inner face of the cloud) of 
—2 

U = 10 (Ferland 1981). This single component model is an obvious 
overs implicat ion . (see Netzer 1982, and Pequignot 1983 ), but the 
results should be of heuristic value. The clouds are ionized by a 
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broken power-law continuum, with a steep continuum (spectra 
index =-1.5) from the optical to 400 ev and a shallower (spectral 
index =-0.5) continuum extending to 100 kev. Metal abundances are 
reduced to 0.5 dex below solar as in Ferland and Netzer (1983)- In 
these models we stop radial integrations when the 
[O I]6300/[0 III]5007 ratio reaches its observed value of 0.05. 

Results of our calculations are shown in Figure 3. The cosmic 


rays 

first 

affect 

the 

emission 

line 

spectrum 

when 

the 

density 

reaches 

f3N* 

-4 

~ 10 

-3 

cm 

As 

was also 

the 

case for 

the 

BLR 

models 

consider 

above 

, the 

first 

lines to be 

affected are 

those produced 

in the 


neutral regions where the photoelectric heating rates are smallest. 

Both the absolute intensity of H«( and the Balmer decrement H*3( /H^*> 

★ 

increase with £N because of the enhanced collisional excitation 

rates caused by the higher temperature. Again, as was true for the 

BLR / lines formed in the ionized regions of the cloud are affected 

relatively less, a consequence of the larger photoelectric heating 

rates there. Eventually the additional heating is large enough to 

keep hydrogen, collisionally ionized and little [O I] emission is ever 
* -1 -3 

formed ( N >10 cm ) - 

Table 2 compares some predictions of models with and without the 

cosmic ray contribution. The cosmic ray density is 

* -i . 5 _3 

£>N =10 cm , the largest value consistent with the presence of 

significant [0 I] emission. Our calculations support the view that 
the relativistic gas associated with the radio lobes could produce 
the high [O III] temperatures in some cases. The ratio 5007/4363 

decreases by a factor of two, indicating a larger electron 
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temperature, when the cosmic rays are added. The temperature of the 

++ 4 * 4 

0 zone doubles from ^1.4x10 K for negligible to 2. 9x10 K 

for £N* = lO -1 * 5 cm -3 - 

The cosmic rays cause a particular enhancement of lines from both 
the neutral and transition regions of the nebula, again because the 
original heating rate was lowest there. In particular, this causes 
the spectrum to resemble that of a LINER in many ways (see Heckman 
1980). A simple analysis of the emission spectrum with the 
additional heat input would deduce sharply enhanced heavy element 
abundances because of the strengthening of the forbidden lines. 


c) Optical Emission from Extended Radio Lobes 

Careful scrutiny of deep optical plates taken of regions within 
the extended radio lobes of AGN often reveals small filaments or 
condensations (cf Graham and Price 1981, Miley et al. 1981, Heckman 
et al . 1982, Burns, Feigelson, and Schreier 1983, Brodie, Konigl, and 
Bowyer 1983). Several mechanisms have been proposed to explain the 
origin of this emission, ranging from shocks to embedded star 
formation, and we wish to consider another. Could these filaments be 
powered by the relativistic gas which surrounds them? 

We model filament 5 of Graham and Price's (1981) paper on Cen A in 
detail. This is a faint condensation roughly 22' away from the 
nucleus along the northeast radio axis, near the center of the middle 
radio lobe . Unlike the cases considered above, there can be little 
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doubt that the optical filaments are in fact surrounded by the radio 
emitting plasma. Other properties, such as the x-ray luminosity, 
distance, and emission measure, are well determined; this object is 
an excellent candidate for detailed study. 

We assume constant gas pressure and an inner gas density of 0.1 
hydrogen atoms per cubic centimeter, resulting in an average density 

_3 

of N "0.2 cm , as deduced by Graham and Price (1981). (Note that 

the cooling time for a gas of this density is /v2000 years, so that a 

heating source is clearly needed to support the observed emission. ) 

For an assumed distance of 3.3 Mpc to Cen A (deVaucouleurs 1979 ), the 

22 

component lies 6.3x10 cm from the nucleus and has a thickness of 
21 

3.2x10 cm. We also assume a uniform heavy element depletion factor 
of 0.1. 


We include the effects of both the ionizing radiation from the 
nucleus as well as the surrounding cosmic rays. A broken power law 
continuum with an x-ray component characterized by a spectral index 
of -0.7, and the steeper index of -1.5 for the optical, joining at 
400 ev, is used. We normalize this continuum to match the 2-6 kev 
flux of 'v 5x10 1( “ ) erg cm ^ s 1 reported by Feigelson et al. (1981). 
The ionization parameter corresponding to this distance, luminosity, 

_3 * ^ 

and density is U = 10 , too small to account for the observed 
ionization and luminosity (as Graham and Price 1981 point out). 

Figure 4 shows the results of varying the cosmic ray density 
between the values where the first effects occur and the upper limit, 
above which no nebular solution is found. Table 3 compares the model 
predictions (both with the ionizing radiation from the nucleus alone, 
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and with both this radiation and a cosmic ray density of 

* -5-65 -3 

= 10 cm ) with the Graham and Price observations. Again, 

the main effect of adding the cosmic rays is to increase the total 
luminosity of the lines by adding more heat and ionization. The 
luminosity increases by an order of magnitude, reflecting the fact 
that the filament is now fully ionized (the model without the cosmic 
rays was only partially ionized) . The relative line strengths do not 
change greatly, although there is a tendency for a higher level of 
ionization because of larger collisional ionization rates. 

The fact that a contribution from; cosmic rays is required to 

reproduce the observed emission line luminosity is the strongest case 

for the fact that the interactions we have considered are actually 

taking place. Although we know of no measurements of the cosmic ray 

density or pressure for the radio iobes surrounding the emission line 

region we consider, we can estimate the cosmic ray density in the 

inner jet (a separation of 4 arc minutes from the nucleus, 
* -5* 9 -3 

N =10 cm ) by combining the estimate of the pressure within 

this jet given by Burns et al. (1983) with the relationship between 
gas pressure and density given in section II. N is probably 
somewhat lower than this in the region we consider, which is nearly a 
factor of 5 more distant from, the nucleus. Considering the wide 
range of possible values for collective heating and , we consider 
this measurement to be in adequate agreement with our model. 
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IV. Discussion 


The calculations presented here have shown that the emission 

line spectra from regions on all distance scales in AGN would be 

affected by the addition of a component of relativistic particles 

with densities similar to those within the radio lobes found on the 

* 

same distance scale. The main effect, for small 0 N , is to increase 

the strength of lines from species with low ionization potentials. 

* —7 

Typically, the effects first become important for N /N e = 10 / 13 and 

* , -4» 4 , 

increase until, for N /N e =10 /A , no nebular solution is found. 
It is important to bear in mind that we have included only heating 
due to direct Coulomb interactions between cosmic rays and cold 
electrons. The scale factor reflects that fact that our results 
are conservative underestimates of the actual effects. Collective 
effects, which increase the heating rate by a factor { 3 , have not 
been included and will enhance the heating rates by many orders of 
magnitude . 

Our original goal in examining this problem was to determine 
whether some properties of AGN's could be explained through the 
addition of a flux of high energy particles. We have discussed 
several aspects of this question above; the causes of some spectral 
differences between radio loud and radio quiet AGN's, production of 
symmetric L<V lines in an outflowing wind, quenching of BLR emission 
in BL Lac objects, the establishment of higher than expected [0 III] 
temperatures in narrow line regions (NLR) of some AGN's, and finally, 
cosmic rays as the energy source for the optical filaments within 
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extended radio lobes. In the last case a conclusive argument for an 
interaction between cosmic rays and the cold gas can be made. 

In principle it should be possible to decide whether the processes 
we have considered are actually taking place in other regions by 
comparing the predicted and observed emission line spectra. 
Unfortunately the differences between nebulae with and without this 
extra heating and ionization are not extreme, and it would be 
difficult to determine whether the additional component of high 
energy particles was present from the emission line spectrum alone. 
For the case of the BLR and NLR models, the addition of the cosmic 
rays actually reduces the agreement between theory and observation 
for some line ratios, while bringing other line ratios into better 
agreement. It is important to remember that our goal here was only 
to determine whether such heat input could be important for several 
standard cases; in particular, we have not reexamined all of 
parameter space to reoptimize model fits. The 'standard' parameters 
we have used in considering- the broad and narrow line regions (such 
as density, composition, ionization parameter) were in turn set by 
optimizing model fits. A different set of parameters would have been 
chosen had the relativistic particles been included in the parameter 
optimization. Nonetheless, it seems difficult to escape the 
conclusion that the additional flux of cosmic rays do little to help 
explain the Grandi and Osterbrock correlation. 

To summarize our results, we would say that the relativistic gas 
which produces the non-thermal radio emission would certainly alter 
the emission line spectrum if the two gasses mix. We consider this. 
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together with 'the models of the extended emission line regions, to be 
the most important results of this paper. Our intention is more in 
keeping with the philosophy of Daltabuit, . MacAlpine, and Cox (1978) 
in their examination of the effects of shocks on a typical BLR cloud; 
these processes could well occur and should be kept in mind. Whether 
the two gases do in fact mix could perhaps best be determined by 
careful high resolution mapping of the emission line regions by the 
Space Telescope, followed by a detailed comparison with VLBI data. 
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Appendix A 

A Crab Nebula Filament 


Henry and MacAlpine (1982) have performed extensive 

predictions of emission line spectra for several filaments in the 

Crab; Nebula. We focus attention on one set of models, those for 

filament 5 of Fesen and Kirshner (1982). As was the case for the 

optical emitters within radio lobes, it is clear that the 

filaments are surrounded by relativistic electrons, although it 

is not clear whether an imbedded magnetic field might shield the 

cld gas from the cosmic rays. Except for the addition of a 

variable component of relativistic electrons, the models we have 

run are identical to model 5.1 of Henry and MacAlpine (1982). 

Basically, it is a constant pressure, plane parallel cloud with 

an inner hydrogen density of 10 ’ cm , a chemical composition 

with enhanced He and C, solar Mg, Si, S, Ar, and Fe, and depleted 

N, 0, and Ne, and ionized 'by a power law continuum extending to 

100 kev. The slab is effectively matter bounded (optically thin 

17 

to hard x-rays), and has a total thickness of 2x10 cm, although 

1 6 

the hydrogen half neutral point occurs at a depth of 6.7x10 cm. 

Tests show good overall agreement between Henry and MacAlpine’ s 

* * 

calculations and ours for small N • 

Figure A1 shows the results of our calculations. We vary the 

* -4 -3 

relativistic particle density between pN =10 cm , where the 

* —1 —3 

first effects occur, and ^N = 10 cm , above which no nebular 

4 6 

solution is found. Two solutions (T V10 K, T/^10 K) can exist 

* -2 -1 -3 

for values of & N in the range 10 - 10 cm , a situation 
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also found in some models of photoionized gas near an x-ray 
source (see Hatchett, Buff, and McCray 1976); we plot the low 
temperature solution. 

The total particle density within the Crab Nebula itself is 

difficult to estimate, but Faraday rotation studies indicate 
-2 -3 

N < 10 cm (Wilson 1974). Direct integration of the spectrum 

8 20 

of the Crab Nebula between 10 and 10 hz indicates 
_6 -5 -3 

N/u lo - 10 cm , depending on the magnetic field. The cosmic 

ray density is probably too small to have any observable effect, 

2 

even if the gasses do mix, unless Q> 10 . 
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Appendix B 

How Reliable are Photoionization Calculations? 

In a series of papers, the most recent of which is Hubbard and 
Puetter (1983), Puetter and co-workers have questioned many of 
the basic premises of the type of calculation we have performed. 
Although this work has resulted in some substantial improvements 
in treatment of some aspects of line transfer, particularly the 
form of the escape probability function used, the series has also 
generated Some- confusion because of the radically different 
nature of the models. In particular, until recently all of their 
calculations assumed constant electron temperature. Their most 
recent calculation performs the thermal balance but treats only 
hydrogen. The fact that their treatment of the line transfer 
process is more realistic than ours, while at the same time they 
neglect the presence of the heavy elements entirely and use an 
oversimplified continuum, raises an important question. Do their 
Calculations disagree with those of Kwan and Krolik (1981, with 
whom we find good agreement) because of sophistication or 
oversimplication? It is this question we address here. 

We have performed a set of calculations in which the 

— 5 

abundances of all heavy elements are reduced by a factor of 10 
We model the two cases considered by Hubbard and Puetter (1983) 
in detail, e.g., a lower density model with typical BLR 
parameters, as well as their high density model. We treat the 
hydrogen line transfer with escape probabilities, and iterate on 
the solution while treating the diffuse fields with the so-called 
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' lambda-iteration 1 procedure. Convergence was fairly rapid, 
typically after only three iterations. Our results, together 
with those of Hubbard and Puetter, are summarized in Table B1 . 
The agreement, particularly for the lower density model, is quite 
good (h/5%). The poorer agreement in the high density model, 
which has a much lower ionization parameter and hence larger 
neutral fraction, is largely due to the fact that Hubbard and 
Puetter neglect secondary ionization by suprathermal electrons in 
neutral regions (see Bergeron and Collin-Souf frin 1973). 

The excellent agreement between our current predictions and 
published thermal equilibrium models was commented on above. The 
fact that we again find good agreement with the Hubbard and 
Puetter calculations (in particular for the model with standard 
BLR parameters) is not surprising. In all these caces, it is 
important to remember the point stressed by Davidson (1973, 
1977); no energy conserving calculation car, go far wrong ♦ This 
is because the strength of the strongest emission lines at a 
given point in a model is simply set by the heating rate at that 
point. This in turn is actually set by the shape of the 
continuum and the recombination coefficient, each fairly simple 
and uncontroversial . In fact, tests show that there has actually 
been fairly little change in the predictions of such models as 
those of Davidson and Netzer (1979) when the same assumptions 
concerning the ionizing continuum and thickness of the cloud are 
made. This is particularly true for strong lines; their 
intensities cannot change greatly in an energy conserving model. 
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Why then do the Hubbard and Puetter calculations disagree with 
those of Kwan and Krolik (as well as our present results for the 
case of nebulae with heavy elements)? Tests show that the 
biggest differences are attributable to differences in the 
continuum shape and their neglect of both heavy element opacities 
and cooling. The first two control the heating rate at great 
depths into the cloud (where the Balmer lines are formed) while 
the latter affects the thermal solution throughout. 

The most important result of this comparison is the fact that 
a computationally expedient technique, the mean escape 
probability formalism, is in excellent agreement with the more 
sophisticated procedures followed by Hubbard and Puetter, when 
the same assumptions concerning ionizing continuum, metalicity, 
and the form of the escape probability function are made. Modern 
techniques for computing emergent spectra from photoionized 
nebulae are both reliable and consistent. 
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Table 1 

BLR Cloud Models 




N*=0 

BN* » 10 4 * 5 cm -3 

L(La) + 


7 .22(4-7) 

9.38(4-7) 

T La 


5.26(4-9) 

7.73(4-6) 

T Ha 


6.19(4-4) 

5620 

La 

1216 

1000 

1000 

Ha 

6563 

481 

248 

Hp 

4861 

106 

73 

He 1 

5876 

11 

16 

He I 

10830 

56 

65 

He II 

4686 

3 

3 

C II 

2326 

41 

127 

C III 

1909 

96 

188 

C III 

977 

48 

153 

C IV 

1549 

608 

1115 

N IV] 

1486 

53 

97 

NV 

1240 

48 

87 

0 III] 

1666 

149 

302 

0 IV] 

1401 

65 

120 

0 VI 

1034 

63 

134 

Mg II 

2798 

215 

96 

Fe II 

RMT40 

51 

2 

Si III 

1895 

25 

550 

Si IV 

1397 

45 

125 

S IV 

1406 

7 

25 


+ erg a”* cm - ^ 
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Table 2 
NLR Models 



N* - 0 

£N* = 10" 1 * 5 cm -3 

L(Hp)* 

0.181 

1.42 

n h ** 

2.24(+22) 

2.53(+22) 


La 

1216 

3431 

23,600 

H8 

4861 

100 

100 

Ha 

6563 

279 

787 

2v 


1052 

12,700 

Hel 

5876 

12 

14 

He I 

10830 

81 

98 

Hell 

4686 

25 

10 

CIl] 

2326 

55 

668 

cm 

977 

25 

125 

CIll] 

1909 

320 

1582 

CIV 

1549 

399 

1204 

TniiI 

6584 

78 

387 

Niv] 

1486 

16 

88 

NV 

1240 

8 

28 

I 01 !, 

6300 

70 

59 

on 

3727 

136 

2036 

on 

7325 

5 

96 

oiii 

5007 

1365 

1172 

oiii 

4363 

24 

65 

[Neill] 3869 

105 

142 

Mgll 

2798 

73 

953 

Fell 

RMT40 

3 

47 

Silll 

1397 

34 

168 

SII 

6731 

165 

898 

SII 

4074 

7 

104 

sin 

1 9532 

189 

115 

ArlV 

1 4740 

6 

10 


*erg cm' 

**cm -2 
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Table 3. 

Extended Emission Regions 


Line 


Observed 

N* - 0 

£N* = 10~ 5 * 65 

F(HB) a 

4861 

- 3 

0.12 

1.3 

La b 

.1216 

— 

4770 

11650 

HI 

2v 

- 

1860 

5440 

Ha 

6563 

330 

299 

433 

HB 

4861 

100 

100 

100 

C II] 

2326 

- 

81 

140 

C III] 

1909 

- 

43 

60 

[N I] 

5199 

- 

11 

11 

[N II] 

6584 

90 

75 

94 

[0 I] 

6300 

- 

36 

29 

[0 II] 

3727 

670 

416 

563 

[0 III] 

5007 

30 

115 

59 

[Ne III] 

3869 

- 

46 

31 

Mg II 

2798 

■ - 

57 

120 

[S II] 

6731 

100 

212 

278 

a Unlts of 

10-17 

erg s"l cm~2 arc sec - ^. 




^Relative to HB . 
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TABLE B1 



Model Comparisons 



Ratio 

P* =■ 0.0207 P = 0.414 


HB Cloudy HB Cloudy 


H a /H 0 

8.5 

8.88 

2.5 

3.23 


5.5 

5.94 

1.1 

1.63 

Bac^^a 

0.35 

0.37 

0.13 

0.27 

V"a 

0.07 

0.073 

0.1 

0.09, 


*dynes/cm^ 
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Figure Captions 

Figure 1. Broad Line Region Clouds . This figure shows 
the effects of various densities of relativistic particles 
on the BLR cloud model of Kwan and Krolik (1981). The scale 
factor P is the ratio of total to Coulomb heating, and can 
be as large as 10^. Intensities are expressed in terras of 
power radiated per square centimeter of cloud. The models 
have all parameters the same as the 'standard' model of Kwan 
and Krolik (1981), with the exception of the cosmic ray 
flux. The column density through the model was set by 
stopping radial integrations when the C 11^ 2326/C IV 1549 
intensity ratio reached its observed value of 0.1. Compact 
radio sources, which occur at roughly the broad line region 

radius, typically have relativistic particle densities of 

1 —3 8 —3 

/viO cm to 10 cm 


Figure 2. Ionization Structures . This figure compares the 
ionization and thermal structure of the Kwan and Krolik 
'standard' model (solid line) with the same model but with 

4.5 _3 

an added flux of cosmic rays of density 10 cm . This 
density is near the logarithmic mean of the values 
determined for compact radio lobes of AGN by Burbidge et al. 
(1974). Both models have column densities set by 
integrating through the cloud until the C IlJ 2326/C IV 1549 
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ratio reaches , the observed value of 0.1. This ratio is 

22 —1 

.reached with a. much lower column density (1.01x10 cm as 
23 -1 

opposed to 10 cm for the Kwan and Krolik model) with the 
cosmic rays because of the additional heating the 
relativistic particles provide. The results of these two 

models are compared in Table 1. 


Figure 3 ^. Narrow Line Region Models . This figure show* the 
results of an added . flux of cosmic rays on a model of a 
typical narrow line region. The models had a variable 
column density which was chosen to maintain a constant 
6300/5007 intensity ratio of 0.05, near the observed value. 


Figure j4. Extended Optical Emission . This figure shows the 
effets of a flux of cosmic rays upon a low density filament 
similar to those typically found in radio lobes of active 
nuclei. Parameters where chosen to be similar to those of 
filament 5 of Cen A as measured by Graham and Price (1981). 
This emission component lies & 22 ' away from the nucleus, 
near the center of the northeast middle radio lobe. This 
model has constant gas pressure, a uniform metal depletion 

-3 

factor of 0.1, an inner density of 0.1 cm , and is ionized 

by both a power law continuum (corresponding to an 

-3 5 

ionization parameter of U=10 * ) and the cosmic rays. 


Thr 
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cloud is only partially ionized in the models with small 

cosmic ray densities, and the level of ionization increases 

* —5 

until the thermal runaway occurs for flN — 10 


Figure Al. Crab Nebula Filaments. This figure shows the 
effects of a flux of relativistic particles on model 5.1 of 
Henry and MacAlpine (1982). Intensities are expressed, in 
total power radiated into 4?T ster radians by a filament 
fully covering the source 


O.f ionizing radiation 



-47- 


Gary J • 

University 


Author's Addresses 


Ferland, Department of Physics and 
of Kentucky, Lexington KY 40506 


Astronomy. 


Richard F. MushotzXy , NASA 
Code 661, Greenbelt MD 20771 


Goddard Space Flight Center, 



1 3N*(cm 3 ) 


Figure 1 




1.0 



(r-r 0 ) 


Figure 2 



F (erg cm 







cn 

CD 



10.0 


0.1 





P (erg 


in 10' 







BIBLIOGRAPHIC DATA SHEET 


1. Report No. 2. Government Accession No. 

TM 86110 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Cosmic Rays and The Emission Line Regions of 
Active Galactic Nuclei 

5. Report Date 

May 1984 

6. Performing Organization Code 
661 “ 

7. Author(s) 

G.J. Ferland (U. KY.) and R.F. Mushotzky 

8. Performing Organization Report No. 

9. Performing Organization Name and Address 

Code 661 

Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 

10. Work Unit No. 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

• | 

14. Sponsoring Agency Code 

15. Supplementary Notes 

Accepted for publication in The Astrophysical Journal 


16. Abstract j 

We discuss the effects that the synchrotron emitting relativistic elec- 
trons, which are known to exist throughout many distance scales in active 
nuclei, could have on the emission line regions which characterize these 
objects. Detailed models of both the inner, dense, broad line region and 
the outer, lower density, narrow line region are presented, together with 
the first models of the optically emitting gas often found within extended 
radio lobes. We show that, in all cases, if the relativistic gas which 
produces the synchrotron radio emission is mixed with the emission line 
region gas then significant changes in the emission line spectrum will 
result. 

The effects of the synchrotron emitting electrons on filaments in the 
Crab Nebula are discussed in an appendix, along with a comparison between 
our calculations, which employ the mean escape probability formalism. 


and recent Hubbard and Puetter models 


17. Key Words (Selected by Author(s)) 

BL Lacertae objects - galaxies: 
nuclei - nebulae: Crab Nebula - 
nebulae: general - radiation 
mechanisms - relativity 

18. Distribution Statement 

f 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

U 

U 

52 



’For sale by the National Technical Information Service, Springfield, Virginia 22151. GSFC 25-44 (10/77) 























